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Abstract 

Context. Asymptotic Giant Branch (AGB) stars are typified by strong dust-driven, molecular outflows. For long, it was believed 
that the molecular setup of the circumstellar envelope of AGB stars is primarily determined by the atmospheric C/0 ratio. However, 
recent observations of molecules such as HCN, SiO, and SO reveal gas-phase abundances higher than predicted by thermodynamic 
equilibrium (TE) models. UV-photon initiated dissociation in the outer envelope or non-equilibrium formation by the effect of shocks 
in the inner envelope may be the origin of the anomolous abundances. 

Aims. We aim at detecting (i) a group of 'parent' molecules (CO, SiO, HCN, CS), predicted by the non-equilibrium study of 
Cherchneff (2006) to form with almost constant abundances independent of the C/O ratio and the stellar evolutionary stage on the 
Asymptotic Giant Branch (AGB), and (ii) few molecules, such as SiS and SO, which are sensitive to the O- or C-rich nature of the 
star. 

Methods. Several low and high excitation rotational transitions of key molecules are observed at mm and sub-mm wavelengths with 
JCMT and APEX in four AGB stars: the oxygen-rich Mira WX Psc, the S star W Aql, and the two carbon stars V Cyg and II Lup. A 
critical density analysis is performed to determine the formation region of the high-excitation molecular lines. 

Results. We detect the four 'parent' molecules in all four objects, implying that, indeed, these chemical species form whatever the 
stage of evolution on the AGB. High-excitation lines of SiS are also detected in three stars with APEX, whereas SO is only detected 
in the oxygen-rich star WX Psc. 

Conclusions. This is the first multi-molecular observational proof that periodically shocked layers above the photosphere of AGB 
stars show some chemical homogeneity, whatever the photospheric C/O ratio and stage of evolution of the star. 

Key words. Astrochemistry, Molecular processes. Stars: AGB and post-AGB, (Stars): circumstellar matter, Stars: mass loss, 
Submillimeter 



1. Introduction 

Circumstellar envelopes of Asymptotic Giant Branch stars 
(AGBs) have long been known to be efficient sites of molecule 
formation. While the outer layers of such envelopes experience 
penetration of interstellar UV photons and cosmic rays resulting 
in a fast ion-molecule chemistry, the deepest layers are dom- 
inated by a non-equilibrium chemistry due to the passage of 
shocks generated by stellar pulsation. Dust forms in those inner 
gas layers, still bound to the star, and grains couple to the gas 
to accelerate it, thereby generating stellar wind and mass loss 
phenomena. The described processes greatly modify the abun- 
dances establish ed by the e quilibrium chemistry in the dense, 
hot photosphere (lTsuiilll973h . 



* Postdoctoral Fellow of the Fund for Scientific Research, Flanders 
** Scientific Researcher of the Fund for Scientific Research, Flanders 



For a long time, the gas chemical composition was believed 
to be dominated entirely by the C/O ratio of the photosphere. A 
C/O ratio greater than one implied that all the oxygen was tied 
in CO, leading to an oxygen-free chemistry, whereas a C/O ratio 
less than one meant that no carbon bearing molecules apart from 
CO could ever form in an oxygen-rich (O-rich) environment. 
This picture, based essentially on thermal equilibrium consid- 
erations applied to the gas, has been first disproved by the de- 
tection of SiO at millimeter (mm) w avelength in carbon-rich (C- 
rich) AGBs (Buiarraba l e"taLlll994 . As for O-rich AGBs, CO2 
infrared (IR) transition lines were detected in various objects 
with the Short- Wavelength Spectromet er (SWS) onboard the 
Infr ared Sp ace Observatory (ISO) (e.g.. lJusttanont et anil996l : 
Rv de etal.lll998l) . 

Theoretical modeling describing the chemistry in the in- 
ner wind of the extreme carbon star IRCh-10216 showed that 
the formation of SiO was due primarily to hydroxyl OH re- 
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action with atomic silicon close to the photosphere as a re- 
sult of shock act ivity and therefore non-equilibrium chem- 
istry (IWillacy & C herchneff 1998). Later on, Duari et al. ( 1999) 
showed that CO2 formation in the O-rich Mira IK Tau results 
from the reaction of OH radicals with CO in the shocked regions, 
implying again that non-equilibirum chemistry was paramount 
to the formation of C-bearing species in O-rich Miras. It was 
then recently proposed that the inner wind of AGBs shows 
a striking homogeneity in chemical composition, despite their 
photo spheric C/O r atio and s tage of stellar evolution dCherchnefB 
I2OO6I) . In p articul ar. ICherchn eff (2006) showed that when taking 
shock chemistry into account, molecules such as SiO, HCN and 
CS are present in comparable amount in the inner layers of M, 
S, and C AGBs, whereas specific molecules (e.g. SO and HS for 
O-rich Miras and C2H2 for carbon stars) are typical for O-rich 
or C-rich chemistries. 

In this letter, we present observations carried out with the 
JCMT and the APEX telescope of four AGBs: one O-rich, WX 
Psc, one S star (= C/O « 1), W Aql, and two carbon stars, II 
Lup and V Cyg. We focus on the detection of (sub)mm transi- 
tions of CO, SiO, HCN, CS, SiS and SO in order to confirm or 
disprove the above hypothesis and to check for homogeneity in 
AGB winds. 

2. Observations and line profiles 

The observations were performed in October 2006 with the 15 m 
JCMTQ for V Cyg, WX Psc, and W Aql, and in the period from 
September till October 2006 with the APEXQ 12 m telescope for 
II Lup, WX Psc and W Aql. Due to technical problems with the 
RxB3 JCMT receiver, only low frequency lines within the RxA3 
receiver (211 - 276 GHz) were obtained. For the APEX observa- 
tions, both the APEX-2A receiver (279-381 GHz) and FLASH 
receiver (460^95 GHz and 780-887 GHz) were used. The ob- 
servations were carried out using a position-switching mode. 
The JCMT data reduction was performed with the SPLAT de- 
voted routines of STARLINK, the APEX-data with CLASS. A 
polynomial was fitted to an emision free region of the spectral 
baseline and subtracted. The velocity resolution for the JCMT- 
data equals to 0.0305 MHz, for the APEX to 0.1221 MHz. For 
WX Psc, W Aql, and II Lup the data were rebinned to a resolu- 
tion of 1 km/s, for V Cyg to 0.75 km/s in order to have at least 
40 independent resolution elements per line profile. The antenna 
temperature, T^, was converted to the main-beam temperature 
(Tmb = T\/rjmb), using a main-beam efficiency rjmi, of 0.69 
for the JCMT RxA3 receiver, of 0.73 for the APEX-2A receiver, 
and of 0.60 and 0.43 for the 46 0-495 GHz and 780-887 GHz 
FLASH channels respectively (Guste n et al.ll2006l) . 

The observed molecular emission lines of four AGB stars in 
our sample are displayed in the Figs. [T]-|4l 

In all stars molecular emission lines of CO, SiO, HCN, and 
CS are detected. Thi s confirms the prediction of homogeneity 
bv lCherchnefl (120061) as these species being 'parent' molecules 
which form in the inner layers of the CSE. This can be under- 
stood in terms of the chemistry of these four molecules being 

' The James Clerk Maxwell Telescope (JCMT) is operated by The 
Joint Astronomy Centre on behalf of the Science and Technology 
Facilities Council of the United Kingdom, the Netherlands Organisation 
for Scientific Research, and the National Research Council of Canada. 
Program ID is m06bn03 (34 h). 

^ APEX, the Atacama Pathfinder Experiment, is a collaboration 
between the Max-Planck-Institut fur Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. Program ID 
is ESO.078.D-0534 (44 h). 



determined by shock propagation and not by the photospheric 
C/O ratio and the stellar evolutionary stage. 

SiS is also detected in all stars and we were able to detect 
with APEX the high-excitation SiS(19-18) line in the O-rich WX 
Psc, the C-rich II Lup and the S-rich W Aql. This is in good 
ag reement with recent S iS OSO, JCMT, and APEX observations 
of ISchoier et al.l (l2007h in a large sample of M and C stars, in- 
cluding WX Psc and V Cyg. This implies that SiS also forms 
close to the star, whatever the stage of stellar evolution. 

Both the SO(65 - 64) and the high-excitation SO(10ii - 
lOio) line were detected in O-rich WX Psc. SO was neither 
found in the S-type W Aql, nor in the two carbon stars II Lup and 
V Cyg. SO appears to be typical for O-ric h AGBs only , suppo rt- 
ing the non-detection of SO in C-stars bv lWoods et al.l ( l2003h . 

Both optically thin and optically thick lines occur (e.g., 
i3CO(2-l) versus the ^"^COil-l) line in W Aql, see e.g. Fig.©. 
The line parameters, i.e., the main-beam brightness temperature 
at the line centre (Tmb), the line centre velocity (w*), and half 
the full line width (ve), are obtai ned by fitting the 'sof t parabola' 
line profile function to the data dOlofsson et al.|[T993l) 

r / \ 21 '^z' 

T{v)^T„,h 1^1 , (1) 

where (3 describes the shape of the line. The velocity-integrated 
intensities are obtained by integrating the emission between 
v■^, ± Ve and are listed in Table [T] A value of /3 = 2 repre- 
sents a parabolic line shape, expected for optically thick expand- 
ing spherical envelopes in case the source is much smaller than 
the radio telescope beam size; [3 < results in a profile with 
horns at the extreme velocit y, expected fo r optically thin lines 
when the source is resolved (lMorrislll975h . In our sample, sev- 
eral lines have a /3-value distinctly larger than 2, indicating that 
the line formation occurs partially in the inner region, where the 
stellar wind has not yet reached its full terminal velocity. This 
is notably the case for the CS, SiO, and HCN lines in our sam- 
ple, coiToborating their formation close to the star. The estimated 
terminal velocities for WX Psc, W Aql, II Lup, and V Cyg are 
19.2 km/s, 18.8km/s, 20.1 km/s, and 12.4km/s, respectively. A 
maximum of 5 % difference is found in the derived terminal ve- 
locity values for each target, being within the velocity binsize, 
indicating that all of the detected lines at least partly survive the 
dust formation and wind acceleration processes. An exception is 
SiS(19-18) in WX Psc (16km/s), but particularly the SO(10ii- 
lOio), with an estimated terminal velocity of only 10 km/s, traces 
a much smaller geometrical region. 

3. Excitation analysis 

In case many transitions of an individual molecule can be ob- 
served, it is possible to assess whether collisional or radiative 
excitation mechanisms can produce the observed line intensi- 
ties. While it is well known that CO is formed in both M, C 
and S-type stars and survives dust condensation, it is of interest 
to study the excitation requirements for th e three other 'pare nt' 
molecules SiO, HCN and CS predicted by CherchnefB (l2006h to 
be abundant in the inner winds of M, S and C AGBs. 

For all lines, except for those from the CO molecule, the 
emission distribution is expected to be much smaller than the 
FWHM beam size of the used telescope. To calculate column 
densities we need to correct for the different beam filling factor, 
/ = 6*1/(01 + with 9b and 63 being the FWHM of the 
beam and the source respectively. To do this, we define a beam- 
averaged brightness temperature, Tb, scaled by coiTecting our 
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Figure 1. Molecular emission detected with APEX and JCMT in the O-rich AGB WX Psc. 
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Figures. Molecular emission detected with APEX in the C-rich AGB II Lup. 
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Figure 4. Molecular emission detected with JCMT in the C-rich AGB V Cyg. 
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Table 1. Overview of the velocity-integrated intensities ( J Tmb dv in K km/s) for the observed line transitions. The frequency is 
listed in GHz and the lower energy level in cm~^. The integrated intensity of detected lines with a low S/N -ratio is given between 
brackets. In case of a non-detection, an upper limit on the integrated intensity is computed as 3crx expected linewidth, with a the 
noise on the data. A '— ' indicates lines that were not observed. APEX-data are reported upright, JCMT-data are given in italics. 
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main-beam brightness temperatures for a fictitious 10" FWHM 
source, i.e, Tb^l/fx Tmb- 

Little interferometric data exist for the molecules that we 
have observed in any circumstellar envelope. We note, however, 
that 4"-6" resolution observations of th e HCN J = 1 — line in 
the Mira variables TX Cam and IK Tau ( lMarvell2005h barely re- 
solve the emission distibutions in these objects. Since, first, those 
objects are closer to the Sun than our target stars and, second, our 
higher exciation lines most likely arise from more compact re- 
gions than the 1 — line, the column densities derived are strict 
lower limits. 

Assuming that the lines are optically thin and that the exci- 
tation temperature, T^^x, between upper and lower level is such 
that Tex^Tbg (with T^g the temperature of any background 
source, e.g. 2.7 K), the integration of the standard radiative trans- 
fer equation shows that ( Goldsm ith & Langer..l999.) 



1.67 X lOi" , , 



(2) 



with Nu the column density of the upper transition state in cm^^, 
/i the transition dipole moment in Debye, u the transition fre- 
quency in GHz, Tmb the main-beam temperature in Kelvin and 
V the velocity in km/s. For rotational transitions, the transition 
moment for diatomic and polyatomic molecule is given by 



/X2(J+1,J) 



J+l 

2J + 3 



Me: 



(3) 



where fj,e is the permanent dipole moment of the molecule. 
Using the velocity-integrated intensities of Table [T] we calcu- 
late the resulting upper state column densities (see Table|2|l. 

Assuming purely coUisional excitation, and ignoring radia- 
tion trapping, one c an derive from solving the statistical equilib- 
rium equation that (lTielensll2005l) 



9l_ I "-crit 

9u V nii2 



+ 1 exp{hiy/kTi^iri) , 



(4) 



where Nu, Ni, gu, gi are the upper and lower state column den- 
sities and degeneracies, Tkin the kinetic temperature, ricrit the 
critical density, and riHa '^he hydrogen density. For a multi-level 
system, the critical density is given by (lTielensll2005h 



ul 



(5) 



where A^i represents the Einstein coefficient, and 7^; the coUi- 
sional rate coefficients. Critical densities were calculated at dif- 
fering te mperatures using th e data available in the LAMDA- 
database dSchoieret al.1 '2005) (see Table ID . Applying Eq. dUl 
together with the column densities from Table |2] the minimum 
density requirements for to achieve the observed number- 
density ratios are tabulated in Table[3] Note that the listed values 
are at Tkin = 300 K, and that for lower temperatures the critical 
density increases. 

The gas temperature, density and velocity structure was cal- 
culated in a self-co nsistent way using the GASTRoNOoM-code 
(iDecin et al ] l2006h to determine where in the envelope the gas 
density falls below the density requirements given in Table|3](see 
Fig. Is). The assumed stellar parameters are listed in Table|4] and 
the derived (maximum) radius of the emitting region is listed 
in Table |3] If coUisional excitation is assumed, it appears from 
Table[3]that the 'parent' molecules SiO, HCN and CS are excited 
in the inner (^ 20 R*) and intermediate (^ 70 R* ) regions of the 
circumstellar envelope and trace regions after the dust conden- 
sation zone where they have been injected from deeper layers. 

It is also of interest to consider the case where coUisional 
excitation is ignored, and the molecules are excited by infrared 
radiation from the star. One can derive that (iTielens,.2005i) 



Ni gi ex^{hv / kT^) — 1 



9u 



W 



1. 



(6) 



where W is the geometrical dilution, being (i?»/2i?)^ when 
R^R*. As seen from Table[3] radiative excitation constrains the 
emitting zone for SiO being ^9R*, for HCN < 13 R*, and for 
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Table 2. Velocity-integrated intensities (in Kkm/s), upper state column densities in cm^^ calculated using Eq. ^ and critical 
densites in cm~^ for the detected multiple-transitions 'parent' molecules used in the excitation analysis. In a few cases, extra 
integrated-intensity values (listed in italics) as found in literature are added to perform the excitation analysis. Literature references 
are given in the footnote. 
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5.1 

3.51 X 10^^ 




5.28 X 10^2 


25.7 
2.06 X 10^^ 


31.9'' 
1.51 X 10" 


7.4 
1.41 X 10" 







at40K n~t 3.06 x 10" 5.26 x 10" 8.41 x 10" 1.27 x 10' 5.83 x 10" 1.32 x 10' 1.85 x 10" 3.22 x 10" 5.09 x 10" 
at 100 K Merit 2.00x10" 3.44 x lO" 5.57 x lO** 8.18x10" 4.22 x lO" 9.63 x lO" 1.35 x lO" 2.37 x lO" 3.76 x lO" 
at 300 K ricrit 1.22x10^ 2.10 x 10** 3.44 x 10** 5.12x10^ 2.35x10^ 5.65 x lO*' 8.35x10^ 1.48x10'^ 2.37 x lO'^ 



'"Hujarrabal et al. ( 1994); ISchoier et al. (200^; '^ IWoods et al] ( |2003[) : Bieging et al. ( 2000. 
^ line detected however with small S/N so that the integrated intensity is quite uncertain 



Table 3. For each target, the first row lists the minimum 
number density for riHa in cm^"^ as derived using Eq. ^ at 
Tkin = 300 K and the second row gives the corresponding 
maximum radius for the emission regions obtained using the 
GASTRoNOoM-code. The third row gives the maximum radius 
for the emitting region calculated using Eq. (|6]l. 









SiO 


HCN 


CS 


WXPsc 


R[Eq. 
R[Eq. 


@] 
©] 


5.0 X 10" 
17 R, 
4.5 R. 


4.4 X 10" 
19 R* 
6R, 


9.0 X 10" 
33 R* 
7.5 R, 


WAql 


R[Eq. 
R[Eq. 


©] 
©] 


2.1 X 10" 
40 R. 
9R. 


2.1 X 10" 
40 R, 
10 R. 


6.34 X 10" 
70 R. 
11 R. 


II Lup 


R[Eq. 
R[Eq. 


0] 
©] 


2.7 X 10' 
13 R, 
2.5 R. 


1.3 X 10" 
53 R. 
13 R. 


7.6 X 10" 
70 R* 
10 R* 


VCyg 


R[Eq. 
R[Eq. 


©] 
©] 


2.8 X 10' 
5R. 


8.7 X 10" 
9R. 
4R. 





Table 4. Stellar parameters used as input for the 
GASTRoNOoM-code. The terminal velocity, v^o, is de- 
rived from the CO lines; the stellar radius from the stellar 
luminosity and temperature. The envelope density falls off as 
~ r^^. Literature references are given in the footnote. 



WXPsc WAql II Lup VCyg 



T. [K] 


2000° 


2800 


2400-' 


1900^ 


R. [10" cm] 


5.5 


2.4 


3.8 


5.1 


L. [IO^Lq] 


10" 


6.8* 


8.8^ 


6.3-^ 


[COMH2] [10-*] 


3" 


6" 


8" 


8" 


distance [pc] 


833" 


230'' 


SOO-^' 


310^ 


Rinncr [R*] 


5" 


8"= 


4/ 


2-f 


Woo [km/s] 


18 


17.5 


21 


10.5 


M [lO-"M0yr-^] 




2.5* 


9/ 


1.23 



"[Decin et al.'^^OOT); ^'Ramstedt et al.' ^200^; '"Knapp & Momsl 
0985); '* Bieging et al. (2000); " Danchi et al. ( 1994); ^ Schoier et al.l 

awey, " Scholer & Olofssoa ( 120010 

' refers to the inner region (A)in lDecin et alj ( |2007|) 



CS ^1 1 R*. In general, these regions are smaller than derived in 
case of coUisional excitation, indicating that the radiation field of 
the star does not have enough energy to sustain the excitation. An 



analogous expression as for Eq. (|6]l can be derived for a radiation 
field characterized by a dust temperature Td at the condensation 
radius i?inncr, with the dilution factor then being (i?innei72-R)^. 
Using the dust condensation radii listed in Table |4] and a dust 
temperature of 800 K the derived emitting regions are a factor 
3.2 larger for WX Psc, a factor 4.3 for W Aql, a factor 4 for 
II Lup, and a factor 2 for V Cyg, resulting in similar radii as in 
case of collisional excitation. 

Although the numbers in Table |3]can only be used as rough 
guidelines, they suggest in both cases a sequence in the exci- 
tation pattern, SiO being the species emitting the closest to the 
star, followed by HCN and CS. 

4. Conclusions 

From the above analysis, one can draw the following conclu- 
sions: 

1 . The observations reported in this letter confirm the status of 
'parent' molecules for CO, SiO, HCN, and CS in AGB stars, 
whose observed molecular lines for m close to the star in sup- 
port of the theoretical predictions of ICherchnefj (l2006l) . The 
excitation analysis suggests that SiO is emitted closest to the 
star, followed by HCN and CS. 

2. High-excitation lines of SiS are detected in all stars, imply- 
ing that SiS too forms close to the star, whatever the stage of 
stellar evolution. However, a thorough line analysis is neces- 
sary to prove or disprove that SiS is more abundant in C stars 
than in O-rich Miras, as predicted by Cherchneff (2006). 

3. SO appears to be typical of O-rich AGBs only. With only two 
excitation lines detected, no information can be drawn on the 
locus of SO formation. However, the SO(lOii-lOio) line has 
an estimated velocity of only 10 km/s in WX Psc when the 
wind terminal velocity for this object is 19.2 km/s. This fact 
suggests that the line excitation occurs close to star In any 
case, the inner chemistry of O-rich AGB envelopes appears 
to be as rich, if not richer, than that of C-rich stars. 

A detailed line analysis of the present data coupled to further 
observational campaigns with JCMT and APEX are planned to 
corroborate these results. 
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Figure 5. The structure of the CSE as derived using the GASTRoNOoM-code is shown for the four studied targets WX Psc, W 
Aql, II Lup, and V Cyg. Since the focus is on the inner and intermediate regions, only the first 200 i?* are displayed. Upper panel: 
Estimated temperature profile, middle panel: estimated gas velocity structure, and bottom panel: estimated hydrogen density nn^ . 
The dashed line indicates the du st condensation ra dius, i?inncr- The variable mass-loss rate of WX Psc, mostly noticeable in the 
hydrogen density, is discussed in lDecin et alj (l2007h . 
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